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The Hamilton-Jacobi (HJ) approach for exploring inflationary trajectories is employed in the 
generation of generalised inflationary non-Gaussian signals arising from single field inflation. Scale 
dependent solutions for /nl are determined via the numerical integration of the three-point function 
in the curvature perturbation. This allows the full exploration of single field inflationary dynamics 
in the out-of-slow-roll regime and opens up the possibility of using future observations of non- 
Gaussianity to constraint the inflationary potential using model-independent methods. The distri- 
bution of 'equilateral' /nl arising from single field inflation with both canonical and non-canonical 
kinetic terms are show as an example of the application of this procedure. 



Density perturbations due to primordial, Gaussian 
scalar metric perturbations are a very good fit to the 
observed structure in the Universe. The simplest model 
for seeding the initial, super-horizon curvature fluctua- 
tions required in this picture is one where small perturba- 
tions in a single inflaton field (f> freeze out on large scales 
during the phase of quasi-de Sitter expansion. The per- 
turbations are expected to be very close to Gaussian. 

It is difficult to distinguish between different inflation 
models if only scalar perturbations are measured. How- 
ever if the non-Gaussian contribution from higher order 
effects could be characterised it would allow the breaking 
of the slow-roll degeneracy and resolve detailed informa- 
tion about the shape of the inflaton potential V(<f>) and 
in principle allow the constraining of more complicated 
forms of the action for (j>. 

Over the past decade much work has gone into the 
quantitative prediction of non-Gaussianity in both single 
field pQ , multi-field models (see [2] for a review) and ones 
with non-canonical forms for the action . The formalism 
extends the calculation of the two-point function in the 
gauge invariant curvature perturbation C, to third order. 
The two-point function, in the Fourier expanded pertur- 
bation £(k), defines the second moment of the distribu- 
tion of the perturbations via an isotropic power spectrum 
(C(k 1 )C*(k 2 )) = (2^) 3 (5( 3 )(k 1 + k 2 )P(fci). The three- 
point function vanishes in the Gaussian case and is de- 
fined via an /nl amplitude [1 

(C(k 1 )C(k 2 )C(k 3 )) = \{2^f5^ (ki + k 2 + k 3 ) x (1) 
5 

/NL(fci,&2,fc 3 ) (-Pfei-Pfes +2 perms) . 

Analytical studies of the prediction for /nl have focused 
on particular configurations of the triangle ki + k 2 + k 3 
('equilateral' k\ — fc 2 = fc 3 , 'squeezed' fci <C fc 2 , fe 3 , etc.) 
and on specific models of inflation. The 'squeezed' sig- 
nal from single field inflation obeys a consistency relation 
that means it is always small /nl^C 1, including in the 
non slow-roll limit [3]. This level of non-Gaussianity 
presents an observational challenge for even the most op- 
timistic future survey forecasts. Large levels of /nl are 
instead predicted in e.g. the equilateral limits in single 



field models in the out-of-slow-roll limit or with non- 
canonical kinetic terms where the speed of sound c s for 
the inflaton can be less than unity with respect to the 
speed of light [4 6J, DBI type models [7], and models 
with multiple fields that undergo large accelerations [HJ. 

A number of authors have examined the numerical pre- 
diction in out-of-slow-roll models that result in large /nl 
[TU] . These have focused on ad-hoc modifications of 
the inflaton potential to induce temporary variations in 
the slow-roll parameters e and/or r\ and have only been 
applied to specific cases. The subtleties involved in the 
numerical evaluation of the integrals required in the in-in 
calculation of the three-point function were explored 
in [gHU]. 

This letter introduces a generalised study of single field 
non-Gaussianity employing the HJ trajectory formalism 
This procedure allows for the calculation of /nl 
arising from random trajectories and thus extends the 
HJ approach for constraining the shape of the inflation- 
ary potential to non-Gaussian observations. This type 
of analysis can complement the similar exploration in- 
volving tensor contributions from different inflationary 
trajectories. Two preliminary examples of these explo- 
rations are shown here; the first being an ensemble of 
single field inflation trajectories with a canonical kinetic 
term (c s = 1) and the second being a case for c s =/= I. 
In both cases results for equilateral /nl as a function 
of wavenumber k are presented. A more general study 
and discussion of the method will be reported in [T2]. A 
choice of units such that that = l/87rG = 1 and 
c = 1 is used throughout. 

Hamilton-Jacobi formalism. — The HJ formalism [TTJ 
IT5] allows for the exploration of all possible inflation- 
ary trajectories consistent with a single inflaton evolving 
monotonically in an FRW universe. The natural basis 
for labelling points in the phase-space of possible trajec- 
tories is an infinite series of 'slow-roll' co-ordinates that 
define the distance of any point along the trajectory from 
pure de Sitter evolution. These are defined as 

and ^ = 2 ^7, T777XTV: ( 2 ) 



H{cj)) J '" " H l ' 

where the value of the scalar field (f> is used as an indepen- 
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dent parameter, H is the Hubble rate, and primes denote 
differentiation with respect to (j>. Re-casting the Fried- 
mann equations into a form where all quantities depend 
on <b one obtains 



= -2H'(<f>) , 



[ff'(0)] a -|ff(# 2 = -^), 



(3) 
(4) 



where t is cosmological time and V((j>) is the inflaton 
potential. This allows one to define an infinite hierarchy 
of differential equations whose solutions self-consistently 
determine possible inflationary trajectories without the 
need to specify a potential V(tf>) 
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where I = 0, f,... and rj = A;, £ = A2 and the number 
of e-folds N, defined via dN/dt = H, has replaced t as 
the independent parameter by making use of the relation 
d<t>/dN = -\/2e. The system together with the 

differential equation dH/dN = —He can be integrated 
with random initial conditions to obtain a Monte Carlo 
sampling of inflationary trajectories |13j . The hierarchy 
can be truncated consistently by imposing that A; = 
for all I > ^ max . The resulting solutions are exact but 
only cover a subset of all possible of solutions, limiting 
the structure and complexity of the trajectories. 

Briefly, the scheme imposed here for Monte Carlo sam- 
pling of HJ trajectories is based on a 'hierarchical prior' 
defined as follows; the initial values of the parameters are 
e = 1, and the rest of the slow-roll parameters are ran- 
domly drawn from a uniform distribution with the ranges 
[-0.1, 0.1] for 77, and [-0.1, 0.1] x 10 2 "' for A/, I > 1. The 
condition Aj max = is imposed. The choice for e ensures 
that inflation ends for that trajectory at the point where 
the random boundary conditions are imposed. The sys- 
tem is then integrated back in time for total of number of 
e-folds N that is itself drawn from a uniform distribution 
with range [40, 80]. This completely fixes the inflationary 
model and, by construction, ensures inflation both ends 
and provides the necessary e-foldings compatible with 
observations. Other choices can be made for the priors 
used in drawing random boundary conditions and also in 
the location of the conditions along the trajectory and it 
should be noted that different choices will affect the final 
distributions in the observables. 

Two-point function. — The trajectories define all back- 
ground quantities that are required for computing the 
correlations of £ on super-horizon scales arising from the 
period of inflation. To do this the late-time solution for 
( is obtained by integrating the Mukhanov-Sasaki equa- 
tion [141 115] arising from expanding the action for £ to 
second order in the perturbations (with constant sound 



(3 + e - 2rj) 



rCfc = 0, 



(7) 



speed c s ): 

dN 2 ' y " ' " "'"dN ' a 2 H 2 ' 

where = C(k) and N is again used as the indepen- 
dent variable. Bunch-Davies initial conditions |16j are 
assumed for each at sufficiently early times when 
c s k 3> aH. The initial condition for H is chosen such 
that this condition holds at time N — 0. The system 
is integrated numerically for a range of k values up un- 
til the mode is sufficiently larger than the horizon (i.e. 
c s k <C aH) to ensure convergence of to a constant. 

The quantity P(;(k) — \Ck\ 2 , at a time when c s k <C aH 
is then the power spectrum of primordial, super-horizon 
curvature perturbations that seeds structure formation. 
In the slow-roll regime, where A; < 1, the dimensionless 
power spectrum will approach a scale invariant solution 
fc 3 P(;(k) ~ fc" s_1 . The n s for each trajectory is calcu- 
lated using a second order finite difference scheme around 
a pivot point fc*. The first order slow-roll prediction is 
n s — 1 2n — 4e [T7]. P^(k) can then be compared to ob- 
servations to constrain the possible set of solutions A; (iV) 
or equivalently H(N) and/or V(4>). The latter being the 
fundamental property of interest in determining the na- 
ture of the inflaton. 

Three-point function. — The inflationary trajectories 
can also be used to obtain generalised non-Gaussian pre- 
dictions. The amplitude of the three-point correlation 
function for £ can be obtained by considering the action 
for a perturbed, minimally coupled scalar field up to third 
order in the perturbation [TJ \5\ |B] 
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where overdots denote differentiation with respect to t, 
di = d/dx.i, and d 2 and d~ 2 are the Laplacian opera- 
tor and its inverse. This action differs from ones usually 
presented in the literature in that no field redefinition 
has been introduced in order to avoid assumptions about 
boundary conditions that would lead to unphysical evo- 
lution of some of the terms for generalised trajectories i.e. 
ones where r\y = 2(e — 77) is still evolving when £fc is out- 
side the horizon. This will be elaborated on in [T2]. The 
three-point function at the end of inflation with TV = N e 
can be computed in the interaction picture by integrating 
the interaction Hamiltonian in time [TJ [5] 

fiV e 

-/ I dN({( kl (N e )( k2 (N e )( k3 (N e ),H int (N)}), (9) 
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FIG. 1. A selection of /nl(&) spectra arising from typical 
trajectories for two cases; constant c s = 1 (bottom) and c s = 
l/v5 (top). The /nl in the c s ^ 1 exhibits strong scale 
dependence where as for c s — 1 it is more mild. 



where H mt is the Hamiltonian for the three-point func- 
tion arising from ([8]). 

This calculation can be carried out analytically in the 
slow-roll regime [TJ [5] and in the out of slow-roll picture 
and/or non-canonical field case if the slow- roll parame- 
ters are constant [TS] . When considering the out of 
slow-roll case with general, non-stationary trajectories 
the calculation must necessarily be carried out numeri- 
cally. A numerical calculation of the integrals involved 
in ([9]) has been carried out for specific inflation models 
with features in the potential [§1 QH] ■ Here ^ is inte- 
grated numerically for each Monte Carlo trajectory [T5] 
The equilateral case of k\ = &2 = k^ = k is chosen in 
what follows but the calculation is valid for all configu- 
rations. The aim in this case is to obtain the function 
/ NL (fc). 

From equations ([I]) and this can be re-cast as 
/nl(^) — 2Im(zfe) where the complex mode Zk satisfies 
the following differential equation 



dz k 
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FIG. 2. The values of the tensor-to-scalar ratio r (bottom) 
and /nl (top) plotted against the scalar spectral index n s . 
All quantities are obtained from the numerical spectra at a 
single pivot scale. Some 10 4 trajectories were included in 
the ensemble with c s = 1. The shaded area shows the 95% 
confidence region for n s allowed by the WMAP 9-year results 
[20] , The red (dashed) line shows /nl~ 5(n s — 1)/12 attractor 
expected in the slow-roll regime for equilateral /nl. 
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with initial condition z k — > as N — > — 00 and are 



The numerical integration of the system requires some 
care due to the highly oscillatory nature of the integrand 
at early times and the evolution of slow-roll parameters 
throughout the integration. Convergence of the integrals 
with respect to both early and late time contributions 
has been verified [T2] . 

Two examples of the results obtained from the numer- 
ical procedure discussed above are reported here. The 
first is an ensemble of O(10 5 ) trajectories with canonical 
speed of sound c s = 1. The second is for a case where 
c s 7^ 1. For simplicity a constant value c s = is 
chosen. From a Lagrangian perspective this corresponds 
to changing the coefficient and power of the kinetic term 
to functions of c 2 s . To include c s evolution an extra set 
of slow-roll parameters needs to be considered and con- 
sequently the functions fi are more complicated. 

The integration is carried out with the aid of GPU 
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FIG. 3. /nl plotted against n a (right) and tinl (left) for the 
same 'slow-roll' trajectories as in Fig. [2] but with constant 
c s = l/s/5. The amplitude of the non-Gaussianity is much 
larger than for the c s = 1 case. The complicated shape of the 
attractor is a consequence of the sinusoidal nature of /nl(&) 
which is not well approximated by a power law. 



parallelisation and only requires on the order of a few 
minutes for the /nl at all k samples to be computed for 
the entire ensemble, /nl spectra resulting from typical 
trajectories with n s « 1 for both ensembles are shown in 
Fig. [T] The amplitude of /nl for the canonical case is 
small, in agreement with analytical estimates for single 
field inflation pQ and mild scale dependence is present 
in some of the solutions. The /nl for the non-canonical 
case have a much higher amplitude, in agreement with 
analytical calculations [5] for the constant out-of-slow- 
roll case, and display a sinusoidal scale dependence in 
Ink. 

The phase portrait of the canonical ensemble is shown 
in Fig. [2] where the tensor-to-scalar ratio r and /nl val- 
ues are plotted against the n s for each trajectory. The 
values shown are all obtained from the numerical spectra 
and are sampled at a fixed pivot scale k+ chosen to be 
to the largest mode in our sample corresponding to the 
scale of the horizon today. 

The distribution of r is as expected with a clear slow- 
roll 'attractor', this is due to the fact that for the fc* 
corresponding to large scales today the spectra are sam- 
pled early on in the inflationary trajectories where the 
slow-roll parameters are typically small given our choice 
of hierarchical priors at the end of inflation. A more 
complex picture may arise from a wider prior, increasing 
^max or for a choice of boundary conditions at the start 



of inflation. The results also show a clear attractor in 
the value of /nl(&*) that agrees with the slow-roll limit 
relation /nl^ i|( n s — 1) for the equilateral case pj\ 

The phase picture for the non-canonical case is shown 
in Fig. [3] The amplitude of /nl at the pivot point is 
plotted against n s and the spectral tilt of /nl defined as 
"nl = din |/NL|/dlnfc at the pivot scale. The attrac- 
tor in this case displays a complicated structure. This 
is due to the fact that, for trajectories that result in 
near scale invariant scalar power spectra, /nl is strongly 
scale dependent and is not well described by a power 
law amplitude and index |21j . A better parametrisa- 
tion of /nl in this case may be of the form f^ik) — 
/nl(1 + 9 sin(wlnfc)). Observations from Planck and 
other future surveys should be able to constrain the am- 
plitudes of /nl seen in this case H23 and potentially de- 
tect any scale dependence of the type seen in these solu- 
tions. 

These results show that it is possible to compute /nl 
for thousands of single field inflation models along with 
the usual n s and r. The method agrees with the ex- 
pected results of /nl for canonical single field slow-roll 
inflation, i.e. small and probably unobservable. However, 
allowing non-canonical c s generally results in observable 
levels of /nl with significant scale dependence. Allowing 
for more structure and widening priors used in sampling 
the trajectories may also result in an extended range of 
amplitudes for /nl- The method will be useful for con- 
straining general single field inflation potentials and/or 
Lagrangian densities (see e.g. 23-25]). Further work 
will explore the landscape of models and, in particular, 
extend the approach to different 'shapes' or configura- 
tions for fci, k2, and k% and make detailed comparisons 
to observations. 
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